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ABSTRACT 

Model tests at M.I.T. have indicated that the cam- 
type wave power generator is capable of high efficiencies 
of energy extraction from ocean waves. This report presents 
a design methodology for determining costs of energy produced 
from wave power generators. Feasible designs were developed 
for ocean locations around the coast of the United States 
with the resulting breakeven costs with oil ranging from 
14 to 30 $/BBL. 
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CHAPTER 1 
INTRODUCTION 

Energy from renewable resources offers an alternative to 
conventional fuel-fired energy production systems. Escalating 
fuel prices, political/economic desires for energy indepen- 
dence and reduction in environmental pollution are all 
incentives for developing renewable resources. The cam-type 
wave power generator is one such resource, utilizing the 
energy in ocean waves. 

The idea of extracting the energy in ocean waves by a 
cam-shaped nodding duck was first presented by Stephen Salter 
of the University of Edinburgh in 1974 (ref. 26) . Since 
waves drive each other with high efficiency, Salter believed 
that an optimally shaped device would cause an incoming wave 
to react as though the device was just another wave. 

In deepwater travelling waves each particle moves in a 
nearly circular orbit. At the surface the diameter of the 
orbiting particles equals the wave height. Below the surface 
the orbital diameters fall off exponentially with depth. 

Salter's duck shape is displayed in figure 1-1. The 
wave side of the cam was designed to conform as nearly as 
possible to the diameters of the orbiting particles. The lee 
side of the cam was designed as a cylinder because this shape 
has constant displacement during the rocking of the cam and 
therefore transmits no energy. 
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Figure 1-1 ; Profile of Salter Cam. 
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Salter tested a model of his device and found that it 
was capable of extracting nearly 90% of the wave energy inci- 
dent on its length. Experiments conducted at M.I.T. with a 
similarly shaped cam also yielded high efficiencies. 

Wave energy, like all other renewable resources, is 
characterized by low intensity and random flux. In order to 
predict the useful absorption of a wave energy device, data 
must be collected and analyzed for the desired ocean location. 
The following parameters adequately describe ocean waves for 
purposes of energy collection: distributions and mean values . 

of wave height, wave period, variability of wave period, power 
density and variability of power density. 

Before its economic feasibility may be determined, the 
Wave Power Generator (WPG) must be measured for performance 
in random seas. To accomplish this a model must be developed 
which is capable of dealing with environmental variables like 
wave height, wave period, period variability, power density 
and power density variability, and with design variables like 
cam diameter, energy storage, power rating and separation 
between WPG units. 

Once performance is determined, subsystem components 
for structure, conversion, electrical generation, transmission 
and energy storage may be selected on the basis of feasibility, 
cost and efficiency. 
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The estimates of deliverable power and subsystem costs 
enable determination of the economic viability of the WPG 
with respect to conventional power sources. 

This report deals with the four basic steps outlined 

above : 

(1) Collection and analysis of ocean wave statistics; 

(2) System performance modelling; 

(3) Subsystem component selection; 

(4) Analysis of investment opportunities. 
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CHAPTER 2 

THEORETICAL AND EXPERIMENTAL BACKGROUND 

Experiments on scale models of the cam have been con- 
ducted at M.I.T. and by Stephen Salter of Scotland. The basic 
characteristics of the M.I.T. cam are depicted in figure 2-1. 
This cam was tested at M.I.T. in both a 2-D mode and a 3-D 
mode: in the 2-D mode the cam's length filled the width of a 

narrow tow tank so that waves generated in the tank had a 
crest length equal to the cam's length; in the 3-D mode the 
cam was tested in a larger tow tank allowing separation 
between the ends of the cam and the tank's walls. In both 
cases the cam was rigidly fixed in space, being allowed only 
to roll about the axis of its length. This axis is parallel 
to wave crest. 

The experiments at M.I.T. were conducted by N.B. Davis 
under the direction of Professor A.D. Carmichael and the 
results are documented in a thesis by Davis (ref. 7) . 

The efficiency of the cam as a wave power absorber is 

defined to be the percent of power in a wave of crest length 

equal to cam length which the cam absorbs. Davis tested the 

cam's efficiency as a function of non-dimensional frequency, 

(jj/a/g, for three levels of inertia: low = .060 Kgm , middle = 

2 2 

.063 Kgm and high = .081 Kgm . His resulting efficiency 
curves are graphed in figure 2-2. Also included in this 
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Figure 2-1 ; M.I.T. Cam Geometry 
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figure is the efficiency curve published by Stephen Salter. 

In each of these curves there is a fairly wide band of high 
efficiency levels. Increasing inertia appears to lower the 
frequency of optimal performance. 

When the cam was tested in the 3-D mode, certain fre- 
quency waves produced very large increases in efficiency. 

This efficiency curve along with the associated 2-D curve for 
the cam with middle inertia is plotted in figure 2-3. Besides 
increasing the cam's efficiency, the 3-D mode displayed 

optimal efficiencies at lower frequencies than the 2-D mode. 

* 

Davis defined a non-dimensional cam response parameter, 
6a/H, where 0 = peak-to-peak radian angular displacement, 
a = cam radius and H = wave height. He plotted response as a 
function of frequency for various levels of damping and this 
plot is included in figure 2-4. Superimposed on this plot 
are dashed lines of efficiency which show that the maximum 
efficiency occurs at the damping level which causes the 
response (0a/H) to have a value ~ 1.0. The level of damping 
was found to have no effect on the frequency of optimal 
performance. 

In a paper by K Budal (ref. 2) it is shown that the 
idealized Salter cam construction comprised of a linear row 
of closely spaced cams is capable of absorbing 100% of the 
power incident on its projected length. This helps to 
explain why the 2-D efficiencies in figure 2-2 are so high. 



Salter 
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Figure 2-2: M.I.T. and Salter Optimum Performance Curves. 
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Figure 2-3: Optimal 3-D and 2-D Performance 
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Figure 2-4 ; 3-D cam response and efficiency for 5 damping levels. 
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Budal also describes the possibility that a wave interaction 
effect on wave-absorbers separated by a distance of approx- 
imately one wavelength (which is large compared to the 
absorber's length) will cause dramatic increases in the 
efficiency of each absorber. This may explain why the M.I.T. 
experiments showed an increased efficiency when the cam was 
operating in the 3-D mode. 
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CHAPTER 3 



ENERGY IN WAVES 



3.1 Theoretical Derivation 

The basic requirement of the wave power generator is to 
capture the energy that is transported across the surface of 
the water in the form of waves. Given a particular set of 
wave characteristics the level of power per unit length of 
free surface (kilowatts/meter) that exists in a train of waves 
may be determined. 

Wave power is the product of the energy in a wave times 
the rate at which the energy is being transported. For plane 
progressive waves the average energy per unit length of free 
surface is 



The average rate of energy flux of a plane progressive 
wave across a vertical control surface is 



E = jp gA 



( 1 ) 




( 2 ) 



The rate at which the wave energy propagates is the 



mean energy flux rate (dE/dt) divided by the mean energy per 
unit length (E) ; 
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dE/dt _ 1 w _ ^ ,,, 

E 2 K g ^ ^ 

Equation 3 shows that the energy in a wave propagates 
at the group velocity (C^) . Since the power in a wave is the 
product of energy times rate of energy propagation, 

P = E Cg = jP9a\ 

= (4) 

For all depths of water, 

^g ^ ^p^^ sinH(2Kh) ^ 



But 



r = Si 

p 27T 



( 6 ) 



Therefore, 



Kh 



p = ^ Q Q«^ 9*^ rl ^ . 

16 ^ ^ s 27 t h sinH ( 2Kh) 



1 2 2 1 
32 ^ Pg H,T[i -P 



Kh 



s '■2 sinH(2Kh) 



■] 



(7) 
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1 Kh 

The value of graphed (figure 3-1) 

or tabulated (table 3-1) . 

1 Kh 

Defining P. = [j * sj.nH(2Kh) ' ' 

In British units, 

P ~ .028 T P* KW/Ft (9) 

s 

In metric units, 

p == T P* KW/M 

o 



( 10 ) 










t 




Wave Period (S) 
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TABLE 3-1 

VALUES OF P FOR VARIOUS WATER DEPTHS AND WAVE PERIODS 
* 

Water Depth (M) 





3 


6 


9 


15 


20 


30 


1.5 


.50 


.50 


.50 


.50 


.50 


.50 


3.5 


.60 


.54 


.50 


.50 


.50 


.50 


5.5 


.52 


.59 


.60 


.54 


.53 


.50 


7.5 


.40 


.53 


.55 


.60 


.59 


.54 


S.5 


.33 


.46 


.51 


.56 


.57 


.58 


11.5 


.29 


.40 


.46 


.53 


.54 


.60 


13.5 


.26 


.34 


.41 


.49 


.50 


.56 


15.5 


.22 


.31 


.36 


.44 


.45 


.53 


17.5 


.20 


.28 


.33 


.41 


.43 


.51 


19.5 


.18 


.25 


.31 


.37 


.39 


.49 


21.5 


.16 


.23 


.28 


.33 


.36 


.47 
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3 . 2 Limitations on Wave Power 

Before examining actual wave power statistics from 
various sites it will be worthwhile to survey some of the 
environmental phenomena which limit the development of power. 

Wind speed is the most important factor in the develop- 
ment of waves. Derbyshire , Sverdrup-Munk and Neumann have all 
proposed spectra relating significant wave height to wind 

speed (U) . These spectra show that H is proportional to 

s 

2 2 5 

U or U * (ref. 15). This implies that the highest values 
of wave power will occur in regions or seasons with very high 
continuous wind speeds . 

Fetch, the distance over water that the wind blows in 

an unobstructed manner, is another important element in wave 

power limitation. For a given wind speed significant wave 

height increases as fetch increases. Figure 3-2 displays 

the effect of fetch on H . 

s 

Related to fetch is sheltering. Harbors and other 
sheltered waterways cannot possibly experience the continuously 
high levels of power available in the open ocean. 

In very shallow water locations wave energy may be 
dissipated in wave-breaking and bottom losses. 

This very general survey should aid in predicting 
locations of the highest power densities. Since winds pre- 
vail from the west for the U.S. it is reasonable to expect 
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coastal power densities to be higher on the West Coast than 
on the East Coast, because of the effects of fetch and 
sheltering. Higher power densities are anticipated at deep 
ocean locations rather than coastal locations because of the 
effects of fetch, sheltering and dissipation. Since winds 
are more severe in the winter months higher power densities 
are expected during winter. The actual values of wave power 
at different locations and during different seasons are 
detailed in Chapter 4 and they substantiate these very general 
observations . 

3 . 3 Probability Distributions 

3.3.1 Wave height . In spectral analysis of wave 
records the Rayleigh distribution is generally assumed for 
wave heights (refs. 4, 15) . The simplified form of this 
distribution is 

P(C) = (11) 

where ? = observed wave height/mean wave height. 

A statistical analysis of this distribution was per- 
formed on 1963 samples from statistics in ref. 12. The 
Rayleigh expression was inverted using the computer library 
program MIT-SNAP and a regression line was fitted to the 



) 

\ 



(>■. 



f. . 

r"" 





















’ - :•’ '4 : 



i ■ •;e5fc 




- : - r Jirlw 




27 



2 

inverted distribution. The regression yielded r = .88, with 

2 

small and evenly distributed residuals. Although the r is 
not particularly high, the regression does indicate that the 
Rayleigh assumption is a reasonable approximation to the 
observed data. 



3.3.2 Wave period . If little is known about the wave 
period at a particular location then a good approximation 
(ref. 10) to the distribution can be achieved by using 



= 1 . e-6’5(|)' 
T 



(12) 



This equation yields a distribution which is symmetric about 
its mean. 

The statistics used in estimating power in this report 

all assume that periods are normally distributed with a 

standard deviation A regression line was fitted to the 

inverted normal distribution for the same 1963 samples used 

2 

in the wave height analysis. The value of r was .96 and 
the residuals were small and uniformly distributed. 



3.3.3 Independence of height and period . Wave heights 
and periods are generally assximed to be independent (refs. 20, 
32) . A statistical test for correlation between height and 
period was performed using 1000 observations from joint 
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height-period tables scattered over all seasons at Nags Head, 
North Carolina (data from ref. 32) . Mean wave periods for 
each of ten wave height groups were tabulated and a hypothesis 
test for correlation was conducted using a 5% risk of type I 
error. Since the data was reduced to only ten wave height 
groups the test was performed using manual methods and the 
t=distribution . 



Hq ; p = 0 



: P 7^ 0 



Test for t = r 




1-r 



where r 



2 _ explained variance 
total variance 




n = number of observations 



The assumed form of the regression equation is 



E{T/H) = a + 0H 



S^(T/H) = -(ZT^ - a ZT - 0ZTH) 
n 



2 1 2 - 
= ^(2T^ - TZT) 
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'ZT 
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'n EH 


s 




ZTH 




EHEH^ 


j 








J 



i 


H. 

1 


T. 

1 


HT 

1 


H .T. 
1 1 


TT 

1 


1 


1 


9.6 


1 


9.6 


92.16 


2 


2 


9.2 


4 


18.4 


84.64 


3 


3 


8.5 


9 


25.5 


72.25 


4 


4 


8.2 


16 


32.8 


67.24 


5 


5 


8.3 


25 


41.5 


68.89 


6 


6 


8.8 


36 


52.8 


77.44 


7 


7 


9.8 


49 


68.6 


96.04 


8 


8 


9.0 


64 


72.0 


81.0 


9 


9 


10.6 


81 


95.4 


112.36 


10 


10 


8.4 


100 


84.0 


70.56 


E 


55 


90.4 


385 


500.6 


822.58 



T = 9.04 



a 




90.4' 




'10 55' 




'8.81 




= 












/N 




.500. 6j 




55 385j 




. -041j 



1 



30 



S^(T/H) =^(822.58 - 8.81 (90.4) - .041 (500.6)) 



.563 



sj = ^(822.58 - 9.04 (90.4)) 



536 



2 



1 



563 

536 



05 



r 



r 



224 



t 




10-2 

1-.05 



65 



For a 5% risk of Type I error, t = 2.3. Therefore, accept 
hypothesis ; p = 0, that there is no correlation between 
H and T. 

3.3.4 Wave power . Nath (ref. 20) derived the pro- 
bability distribution for wave power under the assumptions 

2 

that Power (P) ~ H T, which has already been shown to be 

valid, that height can be represented by the Rayleigh distri- 

2 

bution and that T can be represented by the Rayleigh 

, 2 
distribution. This last assumption for T is the same 

distribution presented in section 3.3.2 for use when wave 

period standard deviation is unknown, that is. 
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T 

P(-) = 1 - e 




T 



T 



2 

Let y = H T = power 



y = mean power 



Using the above assumptions, Nath derived the 



following distribution: 




y 



y = standard deviation 



A statistical test of this distribution was performed 
using the same 1963 samples that were used in the height and 
period tests. The distribution was inverted and a regression 



sion strongly suggests validity of this exponential distribution. 

A further test was conducted on standard deviations of 
wave power arrivals. The exponential distribution yields an 
expected value of o^/P = 1.0. Data from Pierson (ref. 22) 
was utilized with 12 monthly summaries at 20 different loca- 
tions. A regression line was fitted on the observed standard 
deviation of power (Cp) vs . the theoretical standard 
deviation (P) . The regression line was 



line was fitted to the inversion with r 



2 



99. The regres- 
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Op = 2.83 + 1.03 P 

2 

with r = .86. 

With values of Up ranging up to 350 KW/M in the data, 
the regression constant (2.83) is negligible. The slope of 
1.03 compares favorably with the expected slope of 1.0. 

These results further validate the use of the exponential 
power arrival distribution. 

It is important to note that the time between observa- 
tions of wave power may play a significant role in the shape 
of the distribution. Most of the records which were used in 
this report are based on monthly siimmaries of observations 
taken every 12 hours. The distribution of power over any 
given 12 hour interval is clearly not as severe as the 
exponential distribution, largely because of the persistent 
nature of weather patterns. 

It is beyond the scope of this report to attempt to 
determine the shape of the power distribution over a 12-hour 
interval, but a confidence interval about the expected value 
of CTp/P can easily be established. Using 150 sets of data 
from the Marine Environmental Service, Canada (ref. 38) , with 
each data set containing 5 observations of wave energy in a 
given 12-hour interval, a mean and standard deviation of 
energy were computed . 
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n = 150 

X = sample mean = a_/E = .156 

£ 

S = sample standard deviation = .132 
a- = — ^ = .0108 

^ /n^ /TT9 

For a 95% confidence interval, 

X ± 1.96 a- = .156 ± 1.96 (.0108) 

= .156 ± .021 

Range: .135 .177 

In other words over any given 12 hour interval the 
standard deviation of wave energy is expected to be 15.6% ± 
2.1% of the mean wave energy. This is obviously much less 
than the expected 100% using the exponential distribution. 
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CHAPTER 4 
WAVE CLIMATOLOGY 

The problem of defining ocean waves statistically is 
not new, yet, it is only within the past decade that acceptable 
techniques for recording data have been developed. The U.S. 
Army Corps of Engineers Civil Engineering Research Center 
(CERC) has compiled the most accurate data available by its 
system of wave gauges around the coast of the U.S. (ref. 32). 
Unfortunately, all of the wave gauges are located within one 
mile of the coast and are, therefore, not very useful in pre- 
dicting open-ocean power densities. The manual Ocean Wave 
Statistics (ref. 12) provides an acceptable compilation of 
joint wave height and period distributions for very large 
open-ocean areas around the world. These two sources and 
references 22 and 38 provided the data for analysis of ocean 
waves around the U.S. 

Since this report deals with a feasibility study, the 
statistics of ocean waves will be most useful if they present 
a very broad picture. Raw data from 17 coastal locations and 
29 open-ocean locations was analyzed and siammarized to present 
as broad a picture as possible, but still to retain enough 
identity to perform calculations necessary for the study. 

Tables 4-1 to 4-5 contain the following variables: 
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Region - North, Middle and South Atlantic and 

North, Middle and South Pacific off the coast 
of the U.S. 

Ocean location - Coastal waters (within 1 mile) 

Open-ocean (about 100 miles off shore) 
Seasons - Dec-Feb; Mar-May; Jun-Aug; Sep-Nov 
Wave Height - H in meters 
Wave Period - T in seconds 
Wave power density - P in KW/M 
Non-dimensional power variability - CTp/P 
Non-dimensional period variability - a^T 

It is worthwhile to point out some of the strengths and 
limitations of these tables. Hundreds of thousands of obser- 
vations from over 300 records were summarized to create these 
tables; therefore, the tables should not be applied to any 
specific location. Data for coastal regions is more accurate 
than data for deep ocean locations. Calculations of power 
density were performed by the methods outlined in Chapter 3. 
The values of wave period (T) in the deep-ocean sections for 
Atlantic regions are probably not very accurate, as these 
observations are based on shipborne wave-recorders. For this 
reason and since wave period remains essentially unchanged 
with depth (ref. 32) only one value of a^/T was calculated 
for each region, based on coastal readings. 
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The graphs in fig. 4-1 show the variation of the mean 
power level over the seasons for the East and West coasts. 



TABLE 4-1 

WAVE CLIMATOLOGY: DECEMBER TO FEBRUARY 





C 

H 


0 A S 
T 


T A 
P 


L 

<Jp/P 


D E 
H 


E P 
T 


0 C 
P 


E A N 
C^p/P 


(?m/T 

T 


North Atlantic 


1.2 


8.5 


8.7 


.92 


; 3.6 


6.7 


54.6 


.90 


.33 


Mid Atlantic 


1.0 


8.3 


5.2 


1.02 


2.9 


6.2 


31.1 


1.02 


.31 


South Atlantic 


.9 


7.7 


4.0 


.80 


2.5 


6.3 


28.0 


.94 


.38 


North Pacific 










5.2 


10.0 


150 


.98 


.30 


Mid Pacific 


1.1 


10.8 


6.1 


.80 


3.3 


10.9 


90 


1.02 


.25 


South Pacific 


.9 


11.9 


6.2 


.70 


2.4 


11.9 


33 


.90 


.25 
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TABLE 4-2 

WAVE CLIMATOLOGY; MARCH TO MAY 





C 

H 


0 A S 
T 


T A 
P 


L 

^p/P 


D_E 

H 


E P 
T 


OCEAN 
. P. <^p/P 




North Atlantic 


1.1 


8.2 


5.6 


.90 


3.1 


6.3 


38.8 


.97 


.28 


Mid Atlantic 


.9 


8.1 


3.3 


.88 


2.6 


6.0 


24.9 1.10 


.27 


South Atlantic 


.8 


7.6 


2.6 


1.10 


2.3 


6.0 


21.5 1.01 


.36 


North Pacific 










2.7 


9.8 


45 


.88 


.34 


Mid Pacific 


1.0 


10.4 


6.9 


.71 


2.4 


9.7 


27 


.91 


.33 


South Pacific 


. 9 


13.3 


6.6 


.65 


2.1 


11.8 


22 


.83 


.29 



TABLE 4-3 

WAVE CLIMATOLOGY: JUNE TO AUGUST 





C 

H 


0 A S 
T 


T A 
P 


L 


D E 
H 


E P 
T 


0 C 
P 


E A N 
<7p/P 


a/T 


North Atlantic 


.8 


7.4 


2.9 


.86 


2.4 


5.7 


21.2 


1.10 


.20 


Mid Atlantic 


.6 


8.0 


1.9 


1.05 


2.2 


5.6 


16.9 


1.15 


.19 


South Atlantic 


.5 


6.7 


1.5 


1.06 


2.0 


5.6 


14.6 


1.06 


.31 


North Pacific 










1.7 


11.2 


17 


1.07 


;29 


Mid Pacific 


.9 


11.6 


6.7 


.70 


1.5 


10.6 


13 


1.09 


.27 


South Pacific 


.8 


14.0 


6.2 


.55 


1.5 


13.5 


12 


.86 


.24 
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TABLE 4-4 

WAVE CLIMATOLOGY: SEPTEMBER TO NOVEMBER 





C 

H 


0 A S 
T 


T A 
P 


L 


D E 
H 


E P 
T 


0 C 
P 


E A N 
^p/P 


cr,p/T 


North Atlantic 


1.1 


8.7 


5.4 


1.0 


3.0 


6.4 


37.1 


.93 


.31 


Mid Atlantic 


.9 


8.5 


3.7 


.97 


2.6 


6.0 


25.0 


1.05 


.27 


South Atlantic 


.8 


7.0 


2.6 


1.15 


2.3 


6.1 


21.6 


.97 


.36 


North Pacific 










3.9 


11.9 


110 


.93 


.25 


Mid Pacific 


.9 


11.8 


7.2 


.55 


3.3 


10.4 


80 


1.04 


. 28 


South Pacific 


.7 


14.4 


5.4 


.60 


2.2 


14.3 


35 


1.15 


.21 



TABLE 4-5 

WAVE CLIMATOLOGY: ALL SEASONS 





C 

H 


0 A S 
T 


T A 
P 


L 

C7p/P 


D E 
H 


E P 
T 


OCEAN 
P ^p/P 




North Atlantic 


1.0 


8.5 


5.2 


.98 


3.2 


6.3 


37.1 


.98 


.31 


Mid Atlantic 


.8 


7.9 


3.1 


1.03 


2.7 


5.9 


25.6 1.08 


.34 


South Atlantic 


.7 


6.7 


2.5 


1.08 


2.4 


6.0 


22.1 


.99 


.34 


North Pacific 










3.4 


11.0 


81 


.97 


.28 


Mid Pacific 


1.0 


10.4 


5.7 


.70 


2.6 


10.3 


52 


1.01 


.29 


South Pacific 


.9 


13.2 


4.9 


.65 


2.1 


13.2 


25 


.94 


.23 
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Deep Ocean 
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Fifflire 4-1 ; Variation over year of mean power non- 
dimensionalized to mean monthly average. 
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CHAPTER 5 

FORMULATION OF MODEL FOR WPG EVALUATION 

In the previous chapters experimental evidence which 
substantiates the viability of the cam-type wave power 
absorber and the environmental phenomena which affect the 
cam's operation were presented. The next step in the feasi- 
bility study is to analyze the cam with all its associated 
parameters to determine how it might operate in a real-world 
environment. This task may be viewed in the pattern of the 
flowchart below: 




Formulation of the Problem 

The cam system may be characterized by length, dia- 
meter, separation between elements, natural frequency of 
oscillation, motion damping, energy storage, a scheme of con- 
version to energy in a useful form and the users ' demand for 
that energy. 
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The cam's environment (the ocean) is essentially 
characterized at any location by distributions of wave 
height, period or frequency of waves, power or energy in 
waves, water depth and distance from land. 

The coupling between the cam and its environment is 
characterized by sets of efficiency curves which are functions 
of cam design and the environment. 

The basic problem is the same for all energy-related 
systems: that is, to optimally design the parameters of the 

cam (at the feasibility-study level of detail) and to locate 
the system so that it delivers the least expensive energy to 
the users. The model must deal with the stochastic nature 
of the environment and be flexible enough to determine the 
effect on performance of variations in system design. 

Collection and processing of real world data and statistics 

Cam efficiency characteristics and environmental 
variables have been presented in previous chapters. 

Formulation of the model 

There is no necessity to deal with both performance and 
costs simultaneously. Performance will be analyzed first. 

The model must be capable of handling discrete-time 
random inputs, a large n\amber of variables and parameters, 
mathematically intractable efficiency distributions. 
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comparisons of alternative systems or operating conditions, 
reproducible experiments, and output statistics in the form 
of summaries with frequency distributions. 

The obvious model type selection for handling the 
above conditions is simulation. The simulation language GPSS 
was chosen for the following reasons: (1) the author's 

familiarity with its use; (2) discrete-time capabilities; 

(3) ability to deal with complicated distributions; (4) many 
commands for assembling statistics; (5) debugging capabilities; 
(6) transaction flow tracking capabilities; and (7) repeat- 
ability of experiments. The model description is included 
in Appendix A and the prograti listing in Appendix B. 

Model Validation 

Several experiments were performed on the model to 
check its logical and statistical validity. From the table 
POAVA the distribution of power available was plotted to 
check its correlation with the input exponential distribution. 
Several efficiency curves were reproduced by simulating non- 
random seas and extracting data from the table EFFIC (effici- 
ency) . Statistics from the table PEROD (period) were 
correlated with several different input period distributions. 
Individual model segments were checked for proper operation 
using the block count table (a standard output with GPSS) to 
verify that transactions (wave power arrivals) proceeded in a 
logical sequence when the model was constrained to force operation 
of each segment. 
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Experiments were designed and performed to test the 
effects of changing environmental conditions and system 
parameters on system output. These experiments with the 
analyses will be discussed in Chapter 6. 
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CHAPTER 6 

SIMULATION RESULTS AND ANALYSIS 

The simulation model was run with a large number of 
different configurations in order to develop design lanes or 
determine the effects on the system of variable cam diameters, 
hydrodynamic efficiency, wave period variability, wave power 
variability, energy storage, power demand, power rating, 
power availability and cam separation. 

6 . 1 Cam diameter - no separation 

Iterations were performed with the model to determine 
the diameter of maximum power output for a range of wave 
period variances. The maximum output diameter is independent 
of wave period variance but is different for each of three 
levels of inertia tested at M.I.T. The results are listed in 
Table 6-1. 

TABLE 6-1 

CHARACTERISTICS OF MAXIMUM POWER OUTPUT DIAMETER 



Inertia 


cj/a/g 


X/D 


T\/r\ 

max 


Low 


. 66 


7.3 


.83 


Middle 


.76 


5.4 


.93 


High 


.64 


7.7 


1.0 



- i 

iuiiil 
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Note that the maximum power output occurs near the 
apex of efficiency curves from the M.I.T. experiments for 
each level of inertia. The high inertia design delivers more 
power with a smaller cam diameter (smaller co /a/g ) and there- 
fore is the optimal choice of the three. The efficiency 
listed in table 6-1 is non-dimensionalized to the efficiency 
of the high inertia value. The actual number value of the 
efficiency is a function of the variability of wave period 
about its mean: this will be discussed in section 6.3. 

The question next arises as to whether or not this 
single cam diameter indeed delivers the most power from the 
waves since wave periods are varying about the mean and optimal 
cam diameter is a function of wave period. The model was 
tested to see if varying the cam diameter about its optimal 
mean would produce any increase in system efficiency. An 
arbitrary entity was devised representing the standard devia- 
tion of cam diameters about the mean value obtained from 
table 6-1. This entity is defined Figure 6-1 shows the 

simulated results of power out and standard deviation of power 
out vs . The results indicate that the maximum power 

output diameter for a system of cams in a given random sea is 
a point design, in other words, all cams in the system should 
have the same diameter, since the mean power decreases and the 
spread of power increases as diameters are increasingly varied 
about the point mean. 
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Mean Power 

- Band of 1 Std. Dev. 




Fi^re 6-1 ; Effect of Varying Cam Diameters about 
Point Maximum. 
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Since it is possible that the minimum cost/KW system 

will have a smaller diameter than the maximum output system, 

off-design performance was analyzed with the model. The 

performance is depicted in figure 6-2, where power (P) is 

non-dimensionalized with the maximum power output / 

and diameter (D) is non-dimensionalized with the diameter at 

which maximum power is delivered (D ) . 

The spread of power delivered changes considerably as 

diameter is reduced from D due to the filtering character- 

max 

istics of the efficiency curve. Because of this, figure 6-2 

must be adjusted to ensure that all points have equal 

availability. The availability, or the percent of time that 

power is equal to or greater than a specific power level, was 

extracted from the simulation model output tables. The 

availability of power at each point of the curve in figure 6-2 

was compared with the availability of the maximum output 

power (P/P„,^ = 1.0). The curve of adjustment factors is 
rria.x 

included in figure 6-3. 

Assuming that the marginal value of increased availa- 
bility is equal to the marginal cost of increased power, the 
adjustment factor curve may be multiplied point by point with 
the off-design power curve to yield the adjusted curve of 
figure 6-4. This adjusted curve is used when cost estimates 
are applied to find the minimum cost diameter. 
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Fi|S:ure 6-2 : Delivered Power vs Cam Diameter - 

No Separation 



Availability (P) /Availability (Pmax) 
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Figure 6-~^ : Availability Adjustment Factor For Off- 

Design Diameters - No Separation 



max 
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Fi^re 6-4 ? Delivered Power Diameter Adjusted 
For Equivalent Availability. 
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6 . 2 Minimum cost/KW diameter selection - no separation 

With the off-design performance of the cam known, all 
that is needed is the corresponding cost curve to find the 
minimum cost/KW design diameter. The basic trade-off involves 
only structural costs and delivered power. According to 
Kyrkos (ref. 17) structural costs are directly related to the 
weight of structure. A useful structural weight estimating 
relationship (from ref. 25) is 

W. = 1.68341 CN^ + 167.1721 CN - 23.283 tons (1) 

where CN = LBZ x 10“^ 

L = length 
B = beam 
Z = depth 

With the Salter cam, structural weight is essentially 

embodied in the circular cylindrical main body. Therefore, 

2 -5 

beam = depth = cam diameter and CN = LD x 10 
For constant length, 

W^ ~ 1.68341 X 10“^° + 1.671721 x lO"^ 






"■ 



'Mi-- 



,r- 



•«:r 






*■ 00 : 









»>) 



.i2t» la 







M « k: 

- - a- 

i 

•3 ML*" S 



astfT 

-il 



a hi 



Q»b 

A-'«i : 



*>D lOl* 



52 



Therefore, off-design cam costs are proportional to D 

The resulting cost curve, non-dimensionalized to the cost of 

the cam at D is drawn in figure 6-5 along with the off- 
ind.x 

design power curve. Dividing the $/$ curve point-by-point 

$/$ 

by the P/P curve yields or non-dimensional $/KW. 

■' max P/P_ 

max 

From the figure, $/KW is minimized at D/D = .60. 

UlaLX 

In other words, the minimum cost/KW system will have a 

diameter which is 60% of the diameter at maximum power output 

where the cam delivers 68% of the power of D at 36% of the 

max 



cost of D, 



max 



6 . 3 Variability of wave period 

The non-dimensional parameter chosen to describe the 
variability of wave period is cr^T. A given cam design will 
deliver less power as the value of cr^T increases, since the 
cam will be caused to operate outside its region of maximum 
efficiency more often with higher a^T. The simulation model 
was run with various values of and various values of T 
yielding the curves in figures 6-6, 6-7, and 6-8. Figure 6-8 
represents the efficiency curve adjusted for equivalent 
availability. 

Typical values of o^T (from the tables in Chapter 4) 
are displayed in figure 6-8 with the resulting 2-D cam 
efficiencies between 60-70%. 
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Fig:ure 6-5 « Minimtun Cost Diameter Selection - No 
Separation. 
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Pifflire 6-6 : Cam Efficiency -- No Separation, 
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Pi^re 6-7 : Availaliility Adjustment Factor Curve 

For Cam Efficiency C-urve. 



Efficiency 
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Figure 6-8 ; Efficiency Curve Adjusted For Equivalent 
Availability - No Separation. 
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6.4 Cam separation 

Budal (ref. 2) has shown that separation may effectively 
increase the efficiency of a wave power absorber by wave 
interaction. He defines a factor, g, where 

2 N _1 

q = [1 + (N-n)cos(nKd siny) J^(nKd)] (2) 

n=l 

where g = effective efficiency multiplier 
N = number of power absorbers 
K = wave number 
d = separation distance 

Y = angle between cam and normal to wave crest 



Davis (ref. 7) displayed the effect of separation by 
drawing efficiency curves for a 2-D mode and a 3-D mode. In 
the 2-D mode maximum efficiency was about 95%, while in the 
3-D mode maximum efficiency was about 175% (figure 2-3) . 
Since the 3-D mode data was taken in an open tank, the cam 
effectively "saw" images of itself caused by the two tank's 



walls. Davis' data was reexpressed to the form q = 



'3-D 

’2-D 



, and 



plotted against Kd in figure 6-9. The separation distance, 
d, was taken to be twice the distance from the cam's center 
to the tank wall, or the distance from cam to image. 
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Kd 



Figure 6 - 9 : Effect of Separation Using MIT 

Experimental Data. 
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There appears to be a discrepancy between the observed 
"q” and the "q" expected from Budal ' s analysis.^ Assuming 
N = 3 and y = 0, Budal 's equation yields, 

q = [1 + |{2 J^(Kd) + J^(2Kd))]"^ (3) 

This equation is plotted as the solid line in figure 
6-10. For the observed data plotted in the same figure, the 
value of Kd was divided by 3 to achieve the visible correlation 
between the observed data and Budal ' s equation. This is the 
discrepancy referred to previously, and it is also the reason 
that the horizontal axis in that figure is left unmarked. 
Because of the discrepancy, conclusions about the optimal 
value of Kd may not be drawn with certainty. However, it may 
be stated that there exists some value of separation distance 
which will yield appreciable increases in cam efficiency. 

With the 2-D cam, efficiency was shown to decrease 
as the value of cr^T increases. This same effect occurs 
with the 3-D cam. Figure 6-11 shows the cam efficiency as 
a function of cr^/T based on a cam of middle inertia designed 
at the maximum power output, that is, at w/a/g = .63. 



Separation Multiplier (q) 
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0 - MIT Observations 
Budal (N=3) 

Figure 6-10 : Separation Effect Comparison Betv/een 

MIT Observations And Budal Equation. 
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Pi^re 6-11 ; Cam Efficiency Curve With Separation 
At Optimal Separation Distance Adjusted For 
Equivalent Availability. 
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6.5 Minimum cost/KW diameter selection - with separation 
It was shown in section 6.2 that the optimum (minimum 
cost/KW) diameter is not the same as the maximum output 
diameter. An analysis similar to that performed on the 2-D 
cam was performed on the 3-D cam for off-design diameters. 

Davis showed that the maximum power output for the 
3-D cam with middle inertia occurs at a value of co/a/g which 
is lower than that for the 2-D cam: 

. TABLE 6-2 

COMPARISON OF MAXIMUM POWER OUTPUT CHARACTERISTICS FOR 

THE 2-D AND 3-D CAMS 



to/a/g X/D 



2-D 




.76 


5.4 


3-D 




.63 


7.9 


shows 


that for a given 


frequency wave optimal 



cam diameter is smaller for the 3-D cam than the 2-D cam. The 

diameter at maximum power output is 32% smaller for the 3-D 

cam, or, D^^^ (3-D) = .68 D^^^ (2-D). 
max max 

Figure 6-12 shows the off-design 3-D performance; 
figure 6-13 shows the availability adjustment factor; figure 
6-14 shows the adjusted off-design performance curve and 
figure 6-15 shows the minimum $/KW diameter. The analysis is 
similar to that performed in section 6.2. 
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D/D 



max 



Pifflire 6-12 ; Delivered Power vs Cara Diameter - 
vVith Separation. 
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64 




Figti-re 6-13 i Availability Adjustment Factor For Off- 
Design Diameters - With Separation. 







! 



1 

i 






k 



I 

I' 



I 




- ^# 



\ 






i. 



\ 



V 




I 



I 



V. 










max 



65 




Figure 6-14 : Delivered Pov/er Diameter Adjusted 

For Equivalent Availability - With Separation. 



p/p And $/$ 
max max 
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1.25 



D/D 



max 



Pifflire 6-15 : Minimum Cost Diameter Selection 

With Separation, 
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Note that the minimum cost diameter for the 3-D cam is 
about .52 D / while for the 2-D cam it is .60 D . Separa- 

IUd.X ' IU3.X 

tion then will play an important role in reducing costs of the 
system for these three reasons : 

1. Higher efficiency for a given cr^T 

2. Smaller maximum output diameter (D ) 

max 

3. Smaller minimum cost diameter ratio (D/D ) 

max 

These observations may be used to provide a rough 

estimate for the value of separation. At a typical value of 

0 ^T = .30, the factor reduction in diameter of maximum 

power output is . 68 . The factor reduction in the minimum 

cost diameter is .52/. 60 = .867. Using these three factors 

and again assuming that structural costs are proportional to 
2 

D / the structural $/KW of the cam with separation is 

(.68)^ (.867)^ = .25 

1.4 

times the structural $/KW of the cam with separation, assuming 
no additional costs for separation. This remarkable savings 
indicates that further research should concentrate heavily in 
the area of cam separation. 
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6 . 6 Energy storage, power rating and availability 

A very definite trade-off exists betweeci demand 
(assumed to equal power rating for this analysis) and 
availability (% of time that power is equal to or greater 
than demand), for any given level of energy storage. A 
family of curves of demand vs . availability may be drawn 
for any level of energy storage using results from simulation 
model runs. As the example in figure 6-16 displays, there is 
a very sharp "knee" in each curve of demand vs . availability. 
Assuming the marginal value of increased availability equals 
the marginal cost of incremental power, an optimal curve may 
be drawn which is tangent to these "knees". 

The example in figure 6-16 is for an energy storage of 
30 hours, that is, the storage system is capable of delivering 
power at rated capacity for 30 hours. Figure 6-17 displays a 
system of these curves for various levels of energy storage. 

6.7 Optimal level of energy storage 

The curves in figure 6-17 were reexpressed to plot 
demand as a function of energy storage with lines of constant 
availability in figure 6-18. The line of maximum value of 
energy storage represents the highest output per unit of 
storage possible for any availability/demand combination. The 
graph shows that for most levels of availability the maximum 
power output per unit cost of energy storage ranges from 3 to 
7 hours of storage . 
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Pi^re 6-16 ; Sample Storage Curve Construction. 
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Demand Fraction of Power Available (DF) 




FiCTire 6-17: Design Lanes For Storage Cap- 

acity (Hours of Storage Required For 
Given Level of Demand) . 
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Fi|S:ure 6-l8 ; Demand Fraction and Availability 
Tradeoff Lanes. 
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6.8 Variability of wave power (cfp/P) 

The design lanes are based on the exponential power 
function/ whose standard deviation (Up) equals the mean power 
(P) , or Cp/P = 1.0. The tables cited in Chapter 4 on wave 
power availability indicate that cfp/P ~ 1.0 for most deep 
ocean and some coastal locations. However/ there are some 
locations where the value of this parameter may be less than 
1.0. For this reason an adjustment factor needs to be applied 
to the design lanes to account for this change in power 
variability. 

Using the simulation model the adjustment factor (K ) 

s 

was determined and it is graphed in figure. 6-19. The procedure 
for using the curve is as follows : 

1. Determine energy storage requirements from previous 
design lanes; 

2. Enter figure 6-19 with the actual value of Op/P for 

the desired location and retrieve K ; and 

s 

Multiply K by the energy storage from the design 
s 

lanes to find the actual energy storage required. 
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Fig;ure 6-19 i Storage Reduction Factor (K ) 
Variability of Wave Power. 
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CHAPTER 7 

FEASIBILITY TRADEOFF STUDIES 

The objectives of any feasibility study include: 

1. To define a set of feasible alternatives which 
meet initial performance requirements; 

2. To determine cost-effectiveness of the alternatives; 
and 

3. To assume definition of the alternatives to the 
level required for first order cost estimates. 

The feasibility design should be balanced. It is 
inappropriate at the. feasibility level to select specific 
subsystem components unless the user specifies such a selection 
as a design requirement. Subsystems should be identified by 
their general performance characteristics. 

The subsystems analyzed in this chapter include 
structure, conversion, generation, transmission and energy 
storage. The effect of seasons on mean power available and 
a process for power rating selection are also discussed. 

7.1 Seasonality effect on selection of mean power level 

The design lanes in Chapter 6 are all based on a mean 
power level (P) . This mean changes from season to season. 

Using figure 4-1 and assigning P , = 1.0, a typical 

seasonal breakdown may be described as follows : 
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Season 

Summer 

Winter 

Fall/Spring 



season 



.60 P 



annual 



1.40 P 



annual 



1.00 P 



annual 



The question arises concerning selection of the design 
P, that is, should the system be designed for the minimum P 
(summer), the maximum P (winter), or the mean P (fall/spring)? 
Designing for the maximum P yields the most power delivered, 
but designing for the minimum P costs the least. 

The lowest $/KW system depends on the tradeoff between 
delivered power and conversion/generation/transmission (CGT) 
costs. The tradeoff is independent of structural costs 
because maximum power output diameter selection is independent 
of seasons. Energy storage costs are assumed to be constant 
over the seasons. 

The mean power delivered varies with season. If the 
system is designed for the minimum P, then during summer months 
it is rated at 100% of the power available, during spring and 
fall it is rated at 60% of the power available, and during 
winter it is rated at 43% of the power available. The compon- 
ents of this system will have a higher utilization than the 
components of a system designed for the maximum P (winter) . 
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Besides this variation in mean deliverable power, the 
availability of the power delivered changes with season for a 
given design. For example, if the system is designed for the 
minimum P and a value of X is assigned to the availability of 
the rated power during summer, then during winter months the 
availability of the rated power will be greater than X because 
the mean power level in the ocean is higher. 

Assuming that the marginal value of availability 
equals the marginal cost of power, then the equivalent power 
delivered for each of three designs (min P, mean P, and max P) 
may be adjusted for availability by using figure 6-18 and 
holding energy storage constant. 

Table 7-1 lists the power rating, the deliverable 
power, the availability adjustment factor and the equivalent 
power delivered for each of the three designs (min P, mean P, 
and max P) over all seasons, non-dimensionalized to the mean 
P design. 

TABLE 7-1 

RELATIVE POWER DELIVERED FROM SEASONAL DESIGNS 

Equivalent 

Delivered Availability Delivered 



Design 


Power Rating 


Power 


Correction 


Power 


Max P 


1.40 


1.11 


.93 


1.03 


Mean P 


1.00 


1.00 


1.00 


1.00 


Min P 


. 60 


.67 


1.25 


.83 
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If conversion/generation/ transmission (CGT) costs 
represent 0% of the total system cost then the max P design 
is optimal since it delivers the highest equivalent power. 

If CGT costs are 100% of the total system cost then the min P 
design is optimal since it has the least cost. 

Using data in table 7-1, the costs of power for CGT 
costs representing varying percentages of total cost were 
calculated and displayed in table 7-2. The minimum cost/KW 
designs are circled in the table. Since CGT costs are expected 
to lie between 10% and 40% of total system cost the selection 
of the annual mean power level for design considerations is 
the optimal selection. 



TABLE 7-2 

RELATIVE $/KW FOR COMPETING DESIGNS 





CGT 


costs 


as % of total 


cost 


Design 


0 


10 


20 30 


40 


Max P 


*3 


1.02 


1.06 1.09 


1.13 


Mean P 


1.00 


0 






Min P 


1.20 


1.16 


1.11 1.07 


1.02 



50 



1.17 

1.00 
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7.2 Structure 

The most critical structural design element for the 
WPG will be longitudinal bending moments. The Load Line 
Regulations (ref. 4) assert that good design practice requires 
that values of length/depth lie between 10 - 13.5 for ocean- 
going vessels. For vessels operating in the Great Lakes 
regions where wave heights are limited, the Regulations 
require L/D less than 13.5 for vessels less than 325 feet in 
length and L/D less than 19 for vessels with lengths greater 
than 600 feet. One of the reasons for these requirements is 
that insufficient depth lacks the necessary moment of inertia 
to give proper longitudinal stiffness. 

The strength of the cam lies in its circular cylindrical 
backbone. Ocean-going vessels midship cross-sections are more 
nearly approximated by box shapes than circles, implying that 
for similarly sized structures the cam will have a lower mid- 
ship section moment of inertia than the ocean-going vessel. 

With this in mind it would therefore be prudent to restrict 
the cam's maximum value of L/D to 10. If interconnections 
between cams are deemed feasible then the total L/D rather 
than the individual cam's L/D must be less than 10. If cam 
interconnections are flexible then the L/D requirement may be 
relaxed since the flexible link will not carry a load. 
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Structural costs are frequently based on structural 
weight (W) . A set of cost relationships based on Maritime 
Administration data (ref. 17) for ocean-going tanker design 
in 1977 follows: 



Material Cost Ml 
Direct Labor Cost M4 
Indirect Material Cost M5 
Indirect Labor Cost M6 
Material Engineering Cost M7 
Overhead Cost ^ 
Total Cost 



530W 
384W 
31. 8W 
69.12W 
10. 6W 
45.31W 



TC = 1070. 83W 



Structural weight (W) is estimated from the following 
equation : 



W 



YCF*^^ L(1008) 



(1.108-.016 L/D) (1.12-.0163 L/D) 
(35.8-L/D) (14+L/D) 



Where Y = 3.8 - 1.1 (L/D) + .1(L/D)^ 
C = 1.0 + 1.32//iT 
F = 3.8 + 2.1 L/1000 
Z = 5.95 L^D(Cj^+.7) x lO"^ 

(L is in feet) 
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For the cam design, = .785. Holding L/D constant 
at 10, the equation reduces to 



W = 1.075 X 10 



4 (1+1^) (3.8 + *®^(L) 



/L~ 



(L & D in feet) 



This equation was computed for a range of diameters 
from 1-50 meters with the resulting weight/length a function 
of D^: 



Weight (tons) 
Length (m) 



.1865 D 



2 



(D in meters ) 



The proportionality constant, .1865, represents the 

average of the constants calculated for the various diameters. 

The standard deviation of the constants calculated was ±3% of 

the mean, verifying the proportionality of structural weight 
2 

with D for the cam. 

The total cost/length for the cam structure is. 



TC = $1070.83 W 

= $1070.83 (.1865 D^) 

2 

= $199.71 D per meter of length 
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7 . 3 Conversion 

The characteristics of the conversion system are 
governed by the following design requirements; 

(1) the cam's rocking motion must be converted to 
rotational motion at speeds compatible with 
generator requirements; 

(2) damping characterized by variable damping coeffi- 
cients must be utilized to achieve high 
efficiencies ; 

(3) the system must operate with variable speed inputs 
characterized by 6a/H = 1.0 for maximum efficiency 

(4) the system must control power surges; and 

(5) the system must operate in an ocean environment. 

Only one conversion method, hydraulic/pneumatic 

satisfies all these requirements. 

In order to fulfill requirement (1) , variable speed 
radial piston, axial piston or vane pumps may be used. The 
slow speed input implies that piston pumps will be best 
suited to the conversion. 

The pump power takeoff will require a fixed reference 
which may be fulfilled by floating an object behind the cam 
or submerging a ballast weight beneath the cam. In either 
case the fixed reference will be attached to the ends of the 
cam and through seals and bearings the cam will rotate about 



the reference. 
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Hydrostatic drives are best suited to applications 
requiring wide ranges of speeds and torques . The basic 
elements of the hydrostatic drive are a pump, a motor and a 
pump driver. Mechanical energy from the cam driver is con- 
verted to fluid energy by the pump and reconverted to 
mechanical energy by the motor for use by the electrical 
generator. Since it is possible to provide infinitely 
variable ptunp and motor displacements, speed ratios may be 
stepless . 

There are four basic pump/motor configurations: 

(1) Constant capacity pump, constant capacity motor: 

This configuration provides a fixed ratio drive speed 

which is proportional to the displacements of the pump and 
motor. If the two units have equal displacement then the 
drive ratio is one and output speed equals input speed. Speed 
variation may be accomplished by using bypass flow in which 
case power is diverted away from the motor and wasted in the 
form of heat. 

(2) Variable capacity pump, constant capacity motor: 

This configuration provides clutching, constant torque 

and variable speed. With the pump at zero displacement, the 
system is clutched and motor speed is zero. An increase in 
pump displacement increases motor speed, up to the design 
limit, with approximately constant torque. 
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(3) Constant capacity pump, variable capacity motor: 

This system delivers constant power over the design 

speed range. Reduction in motor displacement results in 
increased speed but reduced torque. 

(4) Variable capacity pump, variable capacity motor: 

This arrangement is particularly well-suited for 

designs requiring a wide range of speeds and torques. It 
combines the constant power characteristics of system (3) and 
the constant torque characteristics of system (2) . 

Of these configurations system (4) is best suited to 
the WPG. The variable pump allows selection of the optimum 
damping ratio and the variable motor allows operation over a 
wide range of input speeds without wasting power. Use of this 
system implies coupling with a variable speed generator. 

For the constant speed generator application system (2) 
appears best suited. The variable pump could be controlled 
to select optimum damping. A valving mechanism for wasting 
excess energy, or an energy storage link, would be necessary 
because of the constant motor speed requirement. 

In any configuration gearing is possible between the 
fluid motor and the electrical generator. 

Table 7-3 lists the characteristics of hydrostatic 



conversion . 
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Direct mechanical conversion will not be feasible 
because the cam's velocity reaches zero twice with each 
passing wave. In a purely mechanical system the rocking 
motion could be converted to rotational motion with an 
arrangement of ratchet/pawl and worm/crown gear mechanisms. 
Mechanical, infinitely variable speed transmissions, like the 
Voight unit utilized in wind power, are capable of applying 
damping when coupled with field modulated electrical generators 
But with zero velocity twice each cycle a clutch would be 
needed between the drive and the generator to account for 
inertia effects. Since the transmission must be non-slip in 
order to apply appropriate damping, the clutch requirement 
causes the purely mechanical conversion system to be infeasible 

7.4 Electrical Generation 

Acceptable schemes for generation compatible with power 
grids may be generally classified as constant speed constant 
frequency (CSCF) and variable speed constant frequency (VSCF) . 

7.4.1 CSCF Systems . Synchronous and induction 
generators belong to this category. The constant speed 
requirement dictates the use of energy storage between the 
cam's rocking motion and the generator shaft or a valving 
mechanism which wastes excess energy. Suitable control systems 
are required for sensing and reacting to input speed variations 
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7.4.2 VSCF Systems . There are several types of 
variable speed machines in existence today, all of which are 
more costly than the CSCF machines. 

(1) AC-DC-AC. The AC-DC-AC conversion system is an 

AC generator with DC field excitation. The variable frequency 
DC output is rectified and converted to AC through the use of 
an inverter or a motor/generator set. The system rapidly 
reaches stability with the variable speed input. AC-DC-AC 
may be attractive in deep-ocean installations where DC trans- 
mission to shore is most economical. 

(2) Double Output Induction Generator. This machine 
makes use of an induction generator with rotor energy regen- 
eration. Both rotor and stator supply power to the electrical 
grid. Power from the stator which is generated at grid fre- 
quency is directly fed to the grid. Power from the rotor 
which is generated at slip frequencies is first converted to 
DC, then inverted to AC before being fed to the grid. 
Efficiencies greater than 80% can be achieved with input speed 
variations from 20% to 100% of rated speed. 

(3) Pole-changing Generator. This generator electro- 
magnetically alters the field pole pattern according to shaft 
speed. The electrical output is constant in frequency. There 
are no restrictions on input speeds. 
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(4) Field-modulated Generating Systems. These units 
are generally excited by an AC power source fixing the output 
frequency and modulated by rotational frequency. High shaft 
speeds are required. Output efficiencies are comparable to 
synchronous machines and constant over input speed ratios 
of 5:1. 

Comparisons of off-speed efficiencies for several 
generators are displayed in figure 7-1. 

Table 7-4 compares costs and efficiencies for a number 
of generation units. 

7 . 5 Optimal Power Rating 

It was shown in section 7.1 that the WPG will deliver 
the most power over all seasons if the design is based on the 
mean annual power level. But the selected design power level 
may be more narrowly defined than this. 

Power rating can be optimized. If the power rating 
selected is a small fraction of the mean annual power available 
(P) then the WPG will deliver very little power. If the power 
rating is selected as several times P then the system will 
deliver more power but at a higher cost. For the analysis 
that follows the $/KW of the conversion/generation/transmission 
(CGT) subsystem is assumed to be constant. 
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DC Brushless 
DOIG (VSCF) 

AC Induction (CSCP) 




Fraction of Design Speed 



Figure 7-1 : Off-Design Performance of Several 

Electrical Generators. 
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The fraction of total power delivered up to any given 
level of P may be derived from the simulated zero storage line 
in figure 6-17. Expressing power rating as PR = fraction of 
power available, the resulting curve is plotted as the solid 
line in figure 7-2. 

The two dashed lines in this figure are included for 
comparison of simulated results with the theoretically derived 
exponential power distribution. The upper dashed line repre- 
sents the fraction of all wave arrivals for any given PR from 

“PR 

the equation. Fraction = 1 - e (see Chapter 3) . The lower 
dashed line represents the theoretical fraction of power in 
those waves for any PR. The solid line which also represents 
the fraction of power in waves lies above the lower dashed 
line because of the filtering characteristics of the WPG 
efficiency curve. An empirical formula for the solid line 
was derived: 

Y = K (1 - e”^^) , 

where Y = fraction of total power available 



K 



.90 



rM ■ 



Fraction of Power Available 
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Fraction of Arrivals 
Simulated Fraction Of Power 
Theoretical Fraction Of Power 




Figure 7-2 ; Fraction of Power Available in Y/aves up 
to a Certain Level of Power Rating (PR). 
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Once the power rating is selected the fraction of total 
power captured may be calculated. The following assumptions 
make the calculations more tractable : 

(1) the CGT subsystem delivers all the power in waves 
having power levels less than rated capacity; 

(2) for waves having power levels above the rated 
capacity, the system delivers the rated capacity; 
and 

(3) there is a maximum power level above which the cam 
ceases to operate. 

Assumption (3) requires some clarification. The WPG 
has not yet been tested for operation in severe sea states, 
but it is reasonable to expect that there is some maximum sea 
condition above which the cam cannot operate. In this situa- 
tion the cam may be submerged or allowed to float with power 
generation apparatus secured. This maximum condition may be 

expressed as PR . The fraction of power above PR is lost. 
m3.x rnsx 

Let Y = fraction of total power available which is 
captured by the CGT subsystem 
PR = design power rating 
PRmax “ maximum power rating 

Y = (1) all of the power up to rated capacity + 

+ (2) some fraction of the power between rated 
capacity and maximum capacity + 
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+ (3) none of the power above the maximum capacity. 

— PR 

Term (1) above = K(1 - e ) 

Term (3) =0 

Term (2) was derived in its general form. Referring to 
figure 7-3, term (2) graphically is the ratio of Area A/Area 
(A+B) . In numerical form, term (2) is 

K(PR) [e - e max] 

(1 - e“^^max) (1 + PR ) + PR-PR -(1-e"^^) (1+PR) 

max max 

Therefore, 

v=K ( ^ ,K(PR) [.e~^^-e~^\ax] ^ 

X -CN- \ J- C — PP PP — PP 

(1-e max) (1+ max)+PR-PR ^-(1-e ) (1+PR) 

max 

(1) 

Although the value of PR is not known, PR = 3 is 
^ max max 

a reasonable assumption. Referring to figure 7-2, = 3 

corresponds to 95% of all waves and 84% of all power. In 
other words the assumption of = 3 implies that the 

system will be completely shut down for 5% of the time due to 
excessive sea conditions and 16% of the total power available 
will be lost. 

Equation 1 was calculated for = 3 and plotted in 

figure 7-4. The dashed line in this figure represents 
“PR 

Y = K(l-e ), which corresponds to a system which is shut 
down when power levels exceed PR. 
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Fi^re 7-3 ? Graphical Representation of Fraction 
of Power Between Rated and Maximvun Capacity. 
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Figure 7-4 ? Fraction of Power Delivered vVith Pov;er 
Rating = PR and PR_„„ = 3»0 • 
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Since costs are assumed to be linear with size for the 
CGT subsystem, the slope of the "Y" line in figure 7-4, 
dY/dPR, is actually the marginal cost of the CGT subsystem. 



MP 



cgt 



$/KW cgt 

$/KW system - cgt 



= dY/dPR 



dY ^ + [A] [B] - [C] [D] 

■ [A] 2 

— PP —PR 

where A = (1-e max) (1+PR )+PR-PR - (1-e ) (1+PR) 

1113. X Itl3X 

„ -PR -PR >2 -PR, -PR -PR , 

B = K(e -e max) -2KPRe . (e -e max) 



^ , -PR -PR , 2 

C = KPR(e -e max) 



D = -PRe 



-PR 



Values of dY/dPR (eq.2) have been plotted in figure 7-5. 
Example use of figure 7-5; 

(1) Calculate MF . . 

cgt 

(2) With MP . enter the vertical axis and mark the 

cgt 

corresponding point on the line. 

(3) Read down to horizontal axis for the design power 
rating, which is just PR x P. 
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Piaxire 1 - 5 : Design Lane For Selecting Optimal 

Pov/er Rating And Bandwidth Efficiency Given 
Marginal Cost of CGT ‘ Subsystem, 
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(4) Mark the point on the "Y" line which is vertically 
above or below the point on the line and 

read across to the value of Y. 

Numerical example: 

Design power rating = 1.18 P 
Fraction of power collected (Y) = .77 



The design power rating is optimal for the assumptions 
stated previously. This number multiplied by the cost/KW of 
the CGT subsystem yields the actual cost of the subsystem. 

The fraction of power collected, Y, may be regarded as 
a subsystem efficiency, and it will be referred to in Chapter 
8 as the bandwidth efficiency, 

Efficiencies of specific generation, conversion and 
transmission equipments are assumed to be 1.0 in this analysis 
and they must, therefore, be taken into account in the final 
design process. 

The preceding analysis may be utilized to show that 
variable speed generators should be used in the WPG system. 

A constant speed generator will operate between a cut-in and 
a cut-out wave power level. Below the cut-in level no power 
is delivered. Between the cut-in and cut-out levels the 
generator produces constant power equaling the cut-in power 
level . 
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Numerically, the amount of delivered power for the 
constant speed generator without energy storage is 

T, , T “PR% 

Y - K(1 - e ) = 

/ PR “PR >2 

K (PR) (e- -e max) 

— PP — PP 

(1-e max) (l+PR^^^)+PR-PR^^^-(l-e ) (1+PR) 

ruaX lUaX 

Graphically, the constant speed generator's delivered 

power corresponds to the difference between the "Y" line and 
“PR 

the "K(l-e )" line in figure 7-4. Figure 7-6 compares the 

plots of the variable speed generator (Y) with the constant 

“PR 

speed generator (Y - K(l-e )). Clearly, the constant speed 
generator without energy storage cannot compete with the 
variable speed generator. In addition, the spilled power 
with the constant speed generator must be dissipated as heat 
in the hydraulic conversion system which increases conversion 
costs . 

It was pointed out in section 7.3 that energy storage 
for the constant speed generator must include a hydraulic link. 
But hydraulic energy storage is prohibitive in cost (see table 
7-6) . If the hydraulic conversion link is coupled with 
kinetic energy storage, like a flywheel, then the flywheel 
must physically be located in the cam. The input rotational 
speeds of about 1 RPM from the cam to the hydraulic link would 
have to be increased to about 20,000 RPM for efficient flywheel 
operation. This is infeasible. 
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Variable Sneed Generator 



Constant Speed Generator 




Fi^re 7-6 : Comparison of Power Delivered by Var- 

iable Speed and Constant Speed Generators, 
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With variable speed generators the flywheel, or any 
other energy storage device, could be located ashore and many 
WPG units could feed a single storage device. Flywheel speeds 
could be controlled electrically rather than mechanically. 

For these reasons the variable speed generator is 
preferrable to the constant speed generator. From table 7-4 
the brushless DC generator with processing for inversion to 
grid AC appears to be the most feasible generation subsystem. 



7 . 6 Transmission 

The transmission subsystem depends on system configura- 
tion. For distances less than 20 miles, AC electrical trans- 
mission is the most cost effective method. Beyond 20 miles 
from shore, DC transmission is most cost effective. 

It is not necessary to require transmission of 
electricity. Power generated at the cam may be used in 
hydrolysis of water for production of gaseous hydrogen (GH 2 ) 
or liquid hydrogen (LH 2 ) or even some derivative of hydrogen 
like ammonia. GH 2 niay be feasibly piped over long distances 
and LH 2 may be transported by barge. The hydrogen produced 
may be sold commercially or used to run a fuel cell for 
electrical generation. 

Characteristics of the various transmission schemes 
are displayed in table 7-5. 
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7 ♦ 7 Energy Storage 

Since the power in ocean waves varies so considerably 
with time, some form of energy storage may be required in the 
system. As previously mentioned the CSCF design with hydraulic 
pressure energy storage is infeasible. In VSCF systems energy 
storage may be accomplished by means of flywheels, batteries 
or containerized hydrogen (GH 2 or LH 2 ) . 

Comparisons of various energy storage methods are 
outlined in table 7-6. 

For GH 2 storage additional costs are incurred in 
electrolysis and reconversion to electricity. For LH 2 storage 
additional costs are incurred in electrolysis, liquefaction 
and reconversion. These added costs are outlined in table 7-5. 

From table 7-6 the three most likely forms of energy 
storage are lead-acid batteries, modern flywheels and LH 2 . 

LH 2 is included here because it may become the optimal storage 
medium if the WPG is sited at great distances from shore. 

The technology of lead-acid batteries is well known. 

The only major limitation imposed by the use of these 
batteries is the slow charging rates (in hours for the costs 
associated with batteries in table 7-6) . Lead-acid batteries 
with charging rates on the order of minutes do exist, but the 
costs are several times the cost listed in table 7-6. Slow 
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charge batteries may be feasible if electrical inputs come 
from a large number of sources scattered over a large geo- 
graphic area, which may tend to smooth out the power delivery, 

A considerable amount of interest has been generated 
in the past few years over the use of high energy density 
flywheels. Energy densities have been dramatically increased 
through the use of composites like fused silica and protection 
costs have been reduced because at failure stresses the com- 
posites shatter in millions of small pieces instead of large 
steel chunks like flywheels of the past. Run down times for 
composite flywheels in vacuum chambers have been shown to be 
on the order of months. The limitations on flywheel energy 
storage are: 

(1) very high rotational speeds (20/000-25,000 RPM) ; 

(2) vacu\om chamber sealing; and 

(3) inherent rotational inertia. 

The high rotational speeds can be achieved electrically 
by using a multipole generator to drive a permanent magnet 
motor in the flywheel. The high generator/motor pole ratio 
achieves the desired high rotational speeds. The interface 
is electrical, reducing the vacuum chamber sealing complexity. 
During charge the flywheel acts as a motor, and during dis- 
charge it acts as a generator. The use of the multipole 
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generator is an added cost to cost of flywheel storage. The 
rotational inertia limitation can be eliminated by using 
counter-rotating flywheels or by rigidly securing the unit 
to the earth. 

Taking into consideration all capabilities, limitations, 
costs and efficiencies for the various energy storage schemes, 
the modern flywheel appears to be the most viable system for 
the WPG. 

7 . 8 Other subsystems 

Mooring and direction-keeping subsystems cannot be 
analyzed until free-floating tests are conducted on the cam. 

In Ocean Thermal Energy Conversion (OTEC) systems, mooring 
costs typically run about 10% of structural costs; this cost 
will be assvimed for the WPG. 

There will undoubtedly be a reduction in the effective 
mean power density at a particular site due to the inability 
of the cam to maintain itself parallel to wave crests and its 
inability to maintain a perfect fixed reference, assuming the 
fixed reference is not a structure fixed to the ocean floor. 
These effects, however, will not appreciably distort the 
overall design. If the power density is effectively reduced 
then all system components, except structure, will be down- 
sized accordingly. 
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TABLE 7-3 

CONVERSION SYSTEMS (no storage) 





$/KW 


Efficiency 


Pneumatic 


75 


.55 


Oil hydraulic 


64 


.65 





TABLE 7-4 




ELECTRICAL 


GENERATION 


SYSTEMS 




$/KW 


Efficiency 


Synchronous 


60 


.94 


Induction 


39 


.93 


FMGS 


172 


.80 


Pole changing 


900 


.60 


DC/ AC 


133 


.86 


Multipole 


86 


.95 


Rectifier 


30 


.92 
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TABLE 7-5 

TRANSMISSION (100 MW nominal sizing) 



$/KW 1 mile 


$/KW 100 miles 


Efficiency 


DC electric cable 11 


101 


.95 


GH 2 pipeline 30 


131 


.98 


LH 2 barge 


114 


.95 


Fuel cells 


179 


.70 


Electrolyzer 


89 


.84 


Liquefaction 


234 


.80 



Factor increase in cost* for designs other than 100 MW 
1 MW; 4.3 
50 MV'7; 1.4 



500 MW: 



.7 



. £-qJtc 

■ 4 . 
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Lead-acid battery 
Modern flywheel 
LH2 

Oil hydraulic 
Compressed air 
GH2 



TABLE 7-6 


ENERGY 


STORAGE 


$/KW 


Ib/KWh 


37 


67 


80 


20 


82 


.06 


490 


125 


674 


143 


735 


2.9 



Ft^/KWh 


Efficiency 


2.0 


.60 


.67 


.95 


.01 


.35 


11.1 


.60 


15.1 


.50 


.08 


.40 
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CHAPTER 8 
COST ANALYSIS 



8 . 1 Costs 

Capital costs include costs for structure, conversion, 
generation, transmission and storage. Annual costs may be 
expressed as a % of capital costs. Data was assimilated from 
various sources for ocean-based power generating systems to 
produce the following annual costs; 

% of Capital Cost 



Taxes 2 . 5 
Insurance 2 . 0 
Operating and Maintenance 5 . 5 
Depreciation and Overhead 2.0 
Return 11 . 0 

23.0 



The annual fixed charge rate of 23% must be adjusted, 
or annuitized, for the life cycle of the investment. The 
annuitized fixed charge rate (AFCR) = 

r 

1 _ 
n 



1 



(1+r) 
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where r = fixed charge rate 

n = investment life (years) 

The investment life for the WPG is assumed to be 
20 years which yields APCR = .2337. 

8 . 2 Power Analysis 

The WPG may be analyzed by two methods; (1) fuel 
saver credit - a conventional plant is shut down with the 
resulting savings in fuel; and (2) baseload capacity credit - 
a- conventional plant is physically replaced by the WPG. 

8.2.1 Fuel saver credit analysis . The simplest 
application of the WPG is fuel saving. Power from the WPG 
is fed directly into the utility grid and depending on the 
grid load, conventional fuel-fired plants are backed off with 
a savings in fuel resulting. The value of the fuel savings 
to the utility network requires knowledge of the Incremental 
Heat Rate (IHR; BTU/KWh) . The IHR represents the number of 
BTU's of fuel saved with the reduction of one KWh of conven- 
tional generation. IHR varies constantly. It is currently 
about 9500 BTU/KWh for off-peak power, 11000 BTU/KWh mean 
power and 12300 BTU/KWh for peak power. 
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The value of any fuel-less power generator as a fuel 

saver is 



Fuel Saver Value (<:/KWh) = [(Fuel Cost (<?/10^ BTU) ] 



% 



IHR(BTU/KWh) 



10 



6 



( 1 ) 



Figure 8-1 shows this relationship between fuel costs 
and fuel saver values for various IHR's. 

The busbar cost equation is the basic formula for all 
power cost analyses: 



Power Cost 



/A /T^Tux- [Capital Cost ( $/KW) xAFCR (Fraction/yr ) ]xl00<:/$ 
(Availability (Fraction)] x 8760 hrs/yr 

( 2 ) 



To simplify the use of the equation, the following 
abbreviations will be used: 



FSV = Fuel saver value (<:/KWh) 

BLEC = Baseload energy cost (<?/KWh) 

AV = Availability (fraction) 

CC = Capital cost ($/KW) 

AFCR = Annuitized fixed charge rate (fraction/yr) 



tr ?ao 




- 



- 
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I l_ _J -L 1 1 -L 

5 10 15 20 25 30 35 

Cost of Oil (S/BBL) 



FiCT-re 8-1 : 



System Incremental Heat Rate 
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In fuel saver analysis the following variation of 
the busbar equation applies: 



FSV 



CC (AFCR) 
87.6 (AV) 



( 3 ) 



Once the fuel saver value (FSV) is known, the 
corresponding breakeven fuel cost may be determined from 
figure 8-1. 



8.2.2 Baseload capacity credit . If the availability 
of power from the WPG is sufficiently high then base load 
credit may be applied. The pertinent .equation is, 



BLEC = 



CC (AFCR) 
87.6 (AV) 



(4) 



8 . 3 Availability 

At first glance there is no difference between 
equations (3) and (4) . The actual difference lies in the 
rather complex notion of availability as it applies to the 
WPG. In fuel saver analysis it does not matter when the power 
from the WPG is introduced to the grid: as power arrives from 
the WPG, conventional plants are shut down. For this reason 
the availability of the mean power level in the ocean is 
equal to 1.0 for the fuel saver analysis (neglecting adjustment 
factors) . 







^qiakio' 'H 



Ill 



Therefore, AV„ = 1.0. 
f s 

In the baseload analysis the concept of availability 
for the WPG changes: it is the level of power which the WPG 
can be relied upon to deliver for a certain fraction of the 
year. The level of demand achievable and the energy storage 
required for a given value of availability may be extracted 
from figure 6-18. 

It was mentioned in section 7,1 that the mean power 
level changes with seasons. For this reason baseload credit 
may only be applied to the lowest seasonal value, that is, 
to .60 P, 'the level in summer months. Remaining power may 
be assigned fuel saver credit. The amount applied to baseload 
is .60 (A) (DF) , where A = availability from figure 6-18 and 
DF = demand fraction from figure 6-18. The amount applied to 
fuel savings is (l-AFCR) (1- . 60 (A) (DF) ) , where (1-AFCR) adjusts 
for the difference between fuel savings and baseload value. 

The availability in the baseload analysis then is, = 

.60 (A) (DF)+ (1-AFCR) (1-.60(A) (DF) ) . 

There are two adjustment factors which must be applied 
to both definitions of availability; 

(1) A certain fraction of the year must be set aside 
for equipment maintenance and down-time. Conventional plants 
which are mechanically more complex than the WPG assign 
AV = .75. Ocean Thermal Energy Conversion plants which are 
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less complex mechanically than the WPG assume AV = .95. It 
seems reasonable to set aside a mechanical availability of 
.85 for the WPG. 

(2) In Chapter 7 a method was introduced for optimizing 
power rating (figure 7-5) . The bandwidth efficiency, is 

an adjustment factor which must be applied to availability. 

The two adjustment factors may now be applied to 
and AVj^^ (the bars over and ^j^2. following 

equations imply adjusted versions) . 



AV 



fs 



= .85 n 



band 



(5) 



AVbi = .85 (.60(A) (DF) + (1-AFCR) (1-. 60(A) (DF) ) 

( 6 ) 

The value of (A) in equation (6) for use in figure 
6-18 is .75/ since baseload analysis of the WPG should be 
conducted with respect to conventional plants which assume 
AV = .75. 

Substituting equations (5) and (6) into equations (3) 
and (4) yields the relationships for WPG cost analysis: 



FSV 



CC (AFCR) 
87.6 (^fg) 



BLEC 



CC (AFCR) 
87.6 (^j^i) 



( 7 ) 




,:t <1 
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8 . 4 System Design 

8.4.1 Site selection . The first step in the system 
design process is to select the site for operation of the WPG 
from the tables in Chapter 4. Data should be extracted from 
the "All Seasons" table since the WPG delivers the most power 
when designed to mean annual power levels. 

8.4.2 Diameter selection - without separation effect . 

In Chapter 6 it was shown that D , occurs at wv^a/g = . 64 . 

^ max 

2 

Utilizing the dispersion relation, = .2034 T . The 

minimxim cost diameter was shown to be 60% of D delivering 

max 

68% of the power of D (figure 6-5) . Therefore, design 

max 

2 

diameter » .122 T and the power available = 68% of the 
maximum power. This delivered power may be termed a system 
design efficiency and will be defined Therefore, 

*^design ~ *68. 

2 

Structural cost = $199.71 D . Applying the factor 
increase for mooring cost, structural cost = 1.1(199.71)^ 
<(.122 T^)^ = 3.27 where = power delivered to 

user. 

The cam efficiency is extracted from figure 6-8 using 
the site selection parameter a^T. 
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8.4.3 Diameter selection - with separation effect . 

With separation/ D (3-D) = .68 D (2-D) , and the minimum 

max max 

cost diameter occurs at .52 D , delivering 56% of the power 

max 

at D . Therefore D (3 -d) = .68 (.2034 T^) = .1383 T^, 
max max 

2 2 

and the design diameter = .52 (.1383 T ) = .072 T . 

Structural cost with separation = 1.1 (199.71)^ 

^(.072 T^)^ = 1.139 

The cam's efficiency is extracted from figure 6-11 
using the site selection parameter u^/T. 

8.4.4 Conversion . The hydrostatic drive is the 
selected conversion system with a cost of 64 $/KW and 
efficiency of .65. 



8.4.5 Generation . The DC/ AC generation scheme is the 
selected method of electrical generation with a cost of 133 
$/KW and an efficiency of .86. 

8.4.6 Power rating selection . Power rating (PR) is 
obtained from figure 7-5. The bandwidth of useable power 
(riband^ is also determined from figure 7-5. 
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8.4.7 Transmission . DC Cable transmission is the 
selected method having a cost of 101 $/KW with an 
efficiency of .95. 



8.4.8 Storage . The flywheel storage system is the 
most feasible choice for the WPG. Storage is not necessary 
for the fuel saver analysis, but it may be necessary (or 
desireable) in the baseload analysis. The cost of storage is 
80 $/KWh with an efficiency of .95. The actual storage quantity 
in hours is determined by multiplying the number of hours deter- 
mined from figure 6-18 by the storage reduction factor (K ) in 

s 

figure 6rl9 using the site selection parameter o^/P. The fly- 
wheel input/output device is the multipole m/g set having a 
cost of 86 $/KW with an efficiency of .95. 



8.5 



Design Method for Fuel Saver Analysis 
8.5.1 With separation effect. 



^design 
^conv ~ 



'gen 

’^trans 



= .86 

= .95 



n : from figure 6-11 

'cam ^ 

^overall” ’^design ^cam ^conv ’^gen ’^trans 
= .30 n, 

P 



'cam 



del 



’^overall ^ 




3^T9 

i 

<nf|5 
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structure 
$ 

^ / trans , 'gen 

^CGT = ^n— — + 



1.139 

^del 

$, 



conv 



trans 



trans gen 



trans gen conv 



) (PR) 



= 389 (PR) 



Use figure 7-5 to determine PR and 



MF = 389 (PR) 

^ ^Structure 



perform iteration with values of MF_^„ and PR to make this 



equation valid. 



CC($/KW) = 1«139 T ^ 389(PR) 
^del 



AVfs .85 



FSV = 



87.6 (AV^g) 



8.5.2 Without separation effect . The following 
changes are made to the analysis of section 8.5.1: 



'^design 

'^overall ~ *3^ ’^cam 
ricam * figure 6-8 



8.5.3 Example fuel saver calculation. North Atlantic, 



all seasons, with separation. 



T = 8.5 sec 



P = 37.1 KW/M 



o^T = .31 



From figure 6-11, for a^T = .31, = .83 






^del “ =9.24 KW/M 



^structure = = 64 3 5/KW 



“■cGT = m 



band 



Performing iteration in figure 7-5 yields PR = .68, n 



.64. 



M 
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CC 



s + s 

structure CGT 



643 + 389 (PR) = 908 $/KW 



AVfs .85 



FSV = CC 



87.6 (AV^g) 



FSV = 4.49 <?/KWh 



(.85) (.64) = 

(908) (.2337) 
87.6 (.54) 



.54 



= 4.4 9 <:/KWh 



The corresponding breakeven fuel cost for fuel savings 
= 4.4 9 <?/KWh may be determined from figure 8-1. For IHR = 
11000 BTU/KWh the breakeven fuel cost is 4.05 $/10® BTU, or 
26 $/BBL of oil. 



8 . 6 Design Method for Baseload Analysis 

This method is the same as that for the fuel saver 
except that must be calculated using equation 6. Costs 

must be calculated for zero storage and for the highest 
return storage (figure 6-18) . The lower cost design should 
be selected. 



8.6.1 Example baseload calculation . North Atlantic, 



all seasons, with separation: 




i 



» . 3 
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(1) No storage 
T = 8.5 sec 
p = 37.1 KW/M 
a^/T = .31 
0p/P = .98 
Hq = .25 

P, , = 9.24 KW/M 
ael 

Iband ' 

CC = 908 $/KW 

AVbi = .85 (.60(A) (DF) + (1-AFCR) (1-. 60(A) (DF)) 

from figure 6-18, for A = .75 and zero storage. 



DP * .31 



AV 



bl 



.43 



CC(AFCR) 

87.6 (^bl^ 

= 5.63 <:/KWh 



(2) With storage 

From figure 6-18 the storage quantity yielding 
maximum return is 5.26 hours, at A = .75 and DF = .67. From 
figure 6-19 with o /P = .98, K = .982. Actual storage 

It S 

required = K (5.26) = 
s 



5.16 hours. 
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Storage cost = 80 $/KWh 
'^storage ~ 

Storage I/O cost = 86 $/KW 

Vo “ 

60% of power is diverted to storage: 



^storaae = 

storage (. 95 ) ^(. 95 ) 



AYfai = .85 [.60 (A) (DF) + (1-AECR) (1-.60 (A) (DF) ]=.46 



Hq - .6(.30) ^storage ^1/0 -4 (.30) 



= .23 



Pdel = ^0 ^ = 8.71 KW/M 

CC = 908 i~j) + 426 + 95 = 1514 $/KW 
BLEC = =8.78 C/KWh 

87.6(AVj^j_) 

The cost of baseload energy is lower for the system 



without energy storage than for the system with energy storage 
implying that even the best energy storage scheme, the fly- 
wheel, is still too costly for economical use. It is worth 
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noting, however, that the demand fraction (DF) without energy 
storage is .31, while DF with energy storage ig .67. This 
means that in order to achieve a certain level of baseload 
credit, more than twice as many no-storage WPG's would have 
to be assigned to the baseload as would be required with 
storage-configured WPG's. 

8 . 7 Hypothetical Ocean Energy Farm Cost Calculation 



shipment to shore and conversion to electricity in fuel cells. 
The analysis is conducted for fuel saver with separation 
effect and without energy storage. 



The energy farm electrolyzes water to produce LH 2 for 



Assumed ocean parameters: 



T = 8.5 sec 



P = 100 KW/M 



a^/T = .26 



Gp/P =1.0 



$ 

electrolysis 

$ 

liquefaction 



= 89 $/KW, n = .84 



= 234 $/KW, n = .80 




= 179 $/KW, n = .70 




= 114 $/KW, n = .95 



='li6as 
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% = -13 



Pdel = 



_ 1.139(9-5)^ _ 



structure 



13 



= 714 $/KW 



= 1668 (PR) 



= mi (PR) , yields PR = . 



34 

^band ~ 



CC = 714 + 1668 (. 34) = 128J. $/KW 



^^fs = '^band = 



FSV = 9.14 <?/KWh 



With this high cost occurring in such favorable 
ocean conditions it is unlikely that production of LH 2 will 
be feasible. Costs of electrolysis, liquefaction and fuel 
cells must be drastically reduced and efficiencies improved 
before LH 2 production will be feasible. 
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8 « 8 Value of Energy Storage 

Using the calculations in the previous baseload 
credit analysis, the breakeven cost of energy storage may 
be computed. 



CC + $ 

storage 

CC 



$ 



storage 



.46 

.43 



AVbi (with storage) 
^bl storage) 

(908) - 908 = 63 $/KW 



This is the breakeven cost for 5.16 hours of storage 
plus, the cost of the I/O device. Assuming the I/O device 
represents 20% of total storage cost, the breakeven value of 
one unit of storage is. 



63 (.8) 
5.16 



= 10 



$/KWh 



It is unlikely that any energy storage device will 
achieve this low cost of 10 $/KWh. For this reason energy 
storage is considered infeasible for the WPG. 



8 . 9 Structural Cost Reductions 

The optimization methods for cam diameter in Chapter 
6 were based on the cam being built entirely of steel. Calcu- 
lations of costs in the fuel saver mode showed that structural 
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costs are about 70% of total cost even with the separation 
effect. Structural costs must be reduced. 

It should be possible to build a cam having a small 
cylindrical steel strength member which houses machinery and 
a shroud of foam and plastic which forms the cam's shape to 
the design wavelength. 

Figure 8-2 shows a proposed solution for structural 
cost reduction. 

The size of the cylindrical steel core is governed 

by volume requirements of conversion and generation apparatus. 

Let D = diameter of steel core 

core 

^skin ~ diameter of water surface cylinder 



The volume of the steel core per unit length = 

2 

7T/4 Using the M.I.T. cam design, the volume of the 

shroud encompassing the steel core per unit length = 1.15 '' 

‘'^skin ” ^core* 



Let ^conv ~ conversion system volume density (M /KW) 

3 

Kgen ~ generation system volume density (M /KW) 

2 

K , = $/M of fiberglass coating 

^ s s 

3 

Kpoly = $/M of polyurethane foam filler 

P = power to conversion system (KW/M) 

conv 
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D 



skin 



Polyurethane 
Foam Filler 




Cylindrical Machinery 

Steel Core 



Fi^re 8-2 : Proposed Reduced Cost Structural 

Design. 
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The required core diameter is 



D = [- P (K + n K ) ] 

core IT conv ' conv conv gen^ 



Typical values of the above constants are, 



K = .23 M /KW 

conv 



Vn = ^ /KW 



n = .65 
'conv 



It was shown that the cost of a steel structure is 



199.71 D' 



del 



Therefore , 

$ 



90 P 



conv 



core P 



del 



The cost of the foam filler = 



Vly '’Skin - J °Sore> 



r 
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The cost of the fiberglass surface = 



Vass '1-15 



Typically, 



Vly = 28 $/M- 



Kglass ' 



The total material cost of the shroud = 



22-2 - 22 



Assuming that fabrication costs equal material costs (same 
as for all steel structure) 



shroud 



2'22-2 D3kin ' 22 -h 15.5 



64.4 D , . + 31 D , . - 20 P 

skin skin conv 



and the total cost of the structure is. 



’structure 



^conv ^skin ^1 033^.^ 

^del 



$/KW 



( 9 ) 
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Revised versions of figures 6-5 6-15 may now be 

developed for optimizing structural costs using equation 9. 
Substituting parameters and non-dimensionalizing yields. 



70 (^—) + 64.4 (p5_) 2 + 31 



$ _ 


max 


u 

max 


max 


’max 




165.4 





•Equation 10 is plotted in figure 8-4 for the case 
with separation and in figure 8-3 for the case without 
separation. 

In figure 8-3 the minimiim cost structure without 
separation occurs at D/D = .66 with P/P = .79. 

In figure 8-4 the minimim cost structure with separa- 
tion occurs at D/D _ = .60 with P/P = .73. 

XTioL^c in^^c 

It was shown in section 8.4 that the maximum output 

2 

diameter occurs at .2034 T (without separation) and 

2 

at .1383 T (with separation). 

Total costs may now be expressed in terms of T, with 
the added 10% mooring cost. 



$struct separation) = 

1.1[70 P^^^^+64.4[(. 66) (.2034)]^t'* + 31 ( . 66) { . 2034) T^] 



del 
2 



P + 1-28 + 4.58 T 

conv 

^del 



$/KW 



( 11 ) 
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Figure_8^; Minimum Cost Diameter Selection - No 
Separation. 



p/p And $/$ 
max max 
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D/D 



max 



Pifflxre 8-4 : Minimiim Cost Diameter Selection - 

With Separation. 




'XS'nr ■ 
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^struct separation) = 



1.1[70 P „ + 64.4[(.60) (.1383)]^t'^ + ' 31 ( . 60) ( . 1383) T^] 
conv 



77 P + .49 T + 2.83 T 

conv 

^del 



del 

2 



$/KW 



( 12 ) 



8.9.1 Example fuel saver calculation with steel/foam 
construction . North Atlantic, with separation. 

T = 8.5 sec 
P = 37.1 KW/M 
a^T = .31 

^cam “ (figure 6-11) 

^design ~ 
tIq = .32 

''del " V = 11-9 

'’conv ' '■del = 22.4 KSl/m 

Using equation 12, 

^structure = 277 5/KW 

O Q Q 

MFcgt = ITT yields, PR = .52 

^band 



.60 




3^“i. 
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CC = 377 + 389 (.52) = 579 $/KW 



AVfg .85 



„ _ 579(.2337) 
■ 87.6(.51) 



.51 

3.03 <:/KWh 



This compares with 4.4 9 <:/KWh for the all-steel 
construction . 

Fuel saver analysis calculations were performed 
for all six deep ocean locations around the U.S., for all- 
steel and steel-foam construction, with and without the 
separation effect. The results are displayed' in table 8-1. 



TABLE 8-1 

BREAKEVEN FUEL SAVING VALUES FOR VARIOUS WPG CONFIGURATIONS 









Without Separation 


With Separation 








All-Steel 


Steel-Foam 


All- 


-Steel 


Steel-Foam 




T(sec ) 


P(KW/M) 


% 


1 ?4/KWh 


’'o 


1 ?^/KWh 


^0 


1 (^/KWh 
1 


^0 1 


1 <l>/KMh 
1 


North 

Atlantic 


8.5 


37.1 


.20 


!i1.3 


.25 


5.5 


.25 


1 

1.5 


.32 


3.0 


Mid 

Atlantic 


7.9 


25.6 


.18 


|l3.1 


.2k : 


' 6.2 

j 


.21+ 


' 1.8 


.31 


' 3.1 

1 


South 

Atlantic 


6.7 


22.1 


asi 


1 6.k 


.2k 


1.1+ 


.21+ 


1 3.5 


.31 


2.5 


North 

Pacific 


11.0 


81.0 


.22 

i 


12.9 


1 

.29 


1 

6.1 


.28 ^ 


' 1.5 


.36 

1 


3.1 

1 


Mid 

Pacific 


10.3 


52.0 


.J 

1 


15. 


.28 1 
1 


7.2 


.27 

1 


5.6 

1 


1 

.35 1 


1 

3.5 


South 

Pacific 


13.2 


25.0 


.2h 


64.7 


.32 


27.2 


.32 


19.2 


.ll 


10.5 
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8.10 Analysis of Costs 

From table 8-1 the Middle Atlantic region and the North 
Pacific region yield comparable fuel saver values. The 
penalty for selecting high mean period (T) sites is signifi- 
cant; for a constant fuel saver value, period increases from 
7.9 to 11.0 seconds while power density increases from 22.1 
to 81 KW/M. 

The use of steel-foam instead of all-steel construction 
exhibits reduced costs, but a note of caution must be applied; 
Davis (ref. 7) showed that the cam's natural period of 
oscillation may be effectively raised by increasing the level . 
of inertia of the cam. The steel-foam construction may be 
incapable of supporting the ballast at cam extremities for 
high levels of inertia. 

The separation effect provides the most dramatic 
reduction in costs. For most ocean sites this reduction in 
cost is approximately a factor of two. 

Using table 8-1 and figure 8-1 for IHR = 11000 BTU/KWh, 
the lowest breakeven price of oil is 14 $/BBL, corresponding 
to FSV = 2.5 <:/KWh. It is unlikely that the WPG will ever 
achieve a breakeven cost lower than this. A median FSV = 

5.5 <J:/KWh corresponds to a breakeven oil price of about 
30 $/BBL. 








■3iy>Qf 




V , J 
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The WPG is a resource competing not only with conven- 
tional plants but also with other fuel-less generation 
systems like Wind Power and Ocean Thermal Energy Conversion. 
Costs of these competing resources have been estimated 
(refs. 5, 6, 9 , 13, 16 , 21) to have FSV's ranging from 2.8 
to 3.5 <?/KWh. For the WPG to compete in this range it must 
be configured as follows; 

1. Location in the open-ocean (approximately 100 
miles off-shore) with seas characterized by short 
mean wave periods and/or very high power densities 

2. Separation between elements must be feasible; and 

3. A construction less costly than all-steel must be 
utilized, or a means of reducing the effective 
required diameter must be found. 
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CHAPTER 9 

CONCLUSIONS AM) RECOfflylEHRATIONS 

The follov/ing list siimmarizes the most notable 
conclusions reached in this report: 

1. Ocean Y/aves 

a. V/ave heights and periods are independent; 

b. Wave power distributions may be approximated 
by the exponential distribution; 

c. East Coast waves are characterized by short 
wave periods and low power densities, while V/est Coast 
waves are characterized by long wave periods and high 
power densities; 

d. The mean wave period remains essentially 
unchanged over the seasons; 

e. A typical ocean site displays approximately 
mean annual pov;er levels during the spring and fall, 60^ 
of the mean in summer and 140?^ of the mean in winter. 

2, Modelling 

a. Simulation is a valuable tool for dealing 
with the complex V/PG system; 

b. For a given random ocean site there is one 
optimum cam diameter at a particular level of inertia; 

c. Increasing the level of inertia reduces the 
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required cam diameter; 

d. The most cost-effective diameter is rated at 
68/^ of the mean power available without separation and at 
56?^ with separation for the all-steel stmcture. For the 
steel-foam structure the most cost-effective diameter is 
rated at 79/^ of the mean pov;er available without separa- 
tion and at 735^ with separation; 

e. The parameter may be used in determining 

the cam's hydrodynamic efficiency in a random sea; 

f. The separation effect results in higher effi- 
ciency, smaller maximum output diameter (D ) and 

smaller minimum cost/KW diameter ratio (D/D ); 

• luax 

g. The most cost-effective level of energy storage 
for typical availabilities is 3 - 7 hours. 



3 . Feasibility Studies 

_4 

a. Structural costs are proportional to T ; 

b. Direct mechanical conversion systems are 
infeasible ; 

c. Hydrostatic conversion is well-suited to the 

WPG: 



d. DC electrical generation/transmission yields 
best results; 

e. The bandwidth of useable power may be optimized; 




* 
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f. Flywheel energy storage is the most cost- 
effective storage raedixim. 



4. Costs 

a. The annual fixed charge rate is approximately 

23 /-; 



h. The concept of availability for the WPG is 
complicated by the variability of wave power; 

c. Steel-foam construction produces lov;er costs 
than all-steel construction; 

d. Breakeven costs with oil range from 14 to 
30 $/BBL. 



Recommendations 

Accurate open ocean data needs to be collected, 
since power densities within 1 mile of the coast are too 
low to be economical. Data collection should concentrate 
on finding ocean areas characterized by high power densities 
and lov/ mean annual v/ave periods. 

Model tests of free-floating multi-cam units should 
be conducted to determine (l) the feasibility of utilizing 
the separation effect, (2) mooring forces, (3) direction- 
keeping problems and (4) the limiting increase in natural 
period of oscillation with increases in levels of inertia. 



A 
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Low cost structural concepts may be necessary to 
achieve competitive costs. Such designs should be con- 
ducted with prior knowledge of optimum levels of inertia. 

Inexpensive energy storage will increase the value 
of the WPG, but unit costs must be less than about 10 
$/KWh to warrant storage use. 

The WPG will not likely prove to be the panacea 
of world energy problems. Under the most favorable 
assumptions it is barely competitive on today's energy 
market having a minimum breakeven oil price of 14 S/BBL; 
under unfavorable assumptions it is infeasible with 
respect to today's energy costs. Before the WPG is elim- 
inated from feasible contention as an efficient energy 
producer, further investigation is required. 
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APPETOIX A 
MODEL DESCRIPTION 



FUNCTIONS 

1. NORML: Normal distriLution for wave periods 

(argiiment = RN) . 

2. XPDIS: Exponential distribution for wave power 

(argument = RN) . 

3. CAIIEF: Cam efficiency (argument = ^/D) , 

4. GENEP: Generator efficiency (argument = actual power 

rating seen by generator) . 

5. MULTI: Multiplier for increased power due to cam 

separation (argument = cam separation). 

SAVEVALUES 

1. DEIiND: Demand (K?//) . 

2. POWRT: Design generator power rating (KW), 

3t 4, 5, 7, 11: Not used. 

6. Pump efficiency ifo) , 

8. Mean power level (KW), 

9. Standard deviation of wave period (seconds). 

10. Mean v/ave period (seconds). 

12. CoTinter for assembling power from ten cams. 

13 . Remaining storage capacity (KvV). 
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14. Total storage capacity (ICjV). 

15. Accumijilator of power for delivery to generator. 

16. Counter for power sources. 

VARIABLES 

1. LAMDA: ^/D ( wavelength/diameter) 

2. POV^ER: Power available to each cam = power/length x 

length X IJQLTI. 

3. PO?/LN: Power/length. 

4. PERID: Wave period. 

5. PRATE: Power rating actually seen by generator. 

6. OUTGE: Power outage quantity 

7. PffflPLS: Power available after conversion losses. 

8. CAMLS: Power available after cam losses. 

9. GENLS: Pinal power delivered by system. 

10. SUMPO: Power/length delivered by 10 cams. 

11. POV/: Demand coimter. 

12. TEMPO: Counter. 

13 . AVAIL: Power available in the storage. 

14. KAYD: Kd (v/avenumber x separation distance) . 

MATRIX SAVEVALUES 

1. Array (ixlO) of cam lengths (m). 

2. Array (ixlO) of cam diameters (m). 



TXM 



:10^J 
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3. Array (ixlO) of cam separations (m). 



TABLES 

1 . PEROD : 

2. POAVA: 

3. LI® AD: 

4. EPPIC: 

5 . PULOS : 

6. CALOS: 

7 . GELOS : 

8. OUT AG: 

9. STOAR: 



Period distribution. 

Power in wave distribution. 

^/D distribution. 

Cam efficiency distribution. 

Power leaving the conversion system. 
Power leaving the cam. 

Power out of the system. 

Outage distribution. 

Remaining capacity of storage. 



MODEL CONPIGURATION AI® PLEXIBILITY 



The following model entities may be varied to test 
system response to different input configurations: 
PUNCTIONS 

1. Any empirical cam efficiency distribution may be read 

into the model with CAMEP. 

2. Any empirical generator efficiency curve may be read 

into the model with GENEP. 

3. Purther model tests may yield different efficiency 

multipliers than those currently in the model 
and these may be input using the function MULTI. 




h 
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MTRIX SAVEVALUES 

The model is designed to accept 10 different cam 
configurations for each trial run. The design length of 
the total system was arbitrarily assigned to be 1 ELI with 
each cam allotted l/lO of the total system length. Cam 
diameters may be varied as desired, 

SAVEVAIUES 

1, Demand: Power delivered and storage accumulation 

vary with changes in demand levels. 

2, Power Rating: The ability of the system to deliver 

power depends on the generator power rating. This 
number divided into the actual pov/er available to 
the generator determines the input value for the 
generator efficiency curve, 

3, Conversion Efficiency: This is expressed as a constant 

and accounts for power losses after the cam but 
before the generator, 

4, Mean Power Level: Power levels vary from one location 

to another. This number along with the Function 
XPDIS completely describe the power distribution. 

5, Period: The mean and standard deviation of wave period 

vary from one location to another. These variables 
affect optimum cam diameter, storage and power 
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delivered. 

6. Storage Capacity; The mean power delivered and the 

variability of this power are affected by changes 
in storage capacity. 



TABLES 

All significant model output is in the form of tables. 
Table means, standard deviations and cumulative percen- 
tages are all important identifiers of system output. 



SIMULATE 
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APPEI'IDIX B 
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